Exhibit B 



main 

// This file generates synthetic sinusoids in order to analyze the 
characteristic of 
// a mi rcophone. 

#include <cmath> 

#include <cstdlib> 

#include <cstdio> 

#include <iostream> 

#include <fstream> 

#include <ctime> 

#include <algorithm> 

#include <functional> 

#include <complex> 

#include "audio/wavFile .H" 

#i ncl ude "audi o/Audi oBuf f er . H" 

#include "nsp.h" 

//#define nsp_usesLms 
//#i ncl ude <nsp.h> 



#ifndef M_PI 

#define M_PI 3.14159265358979323846 
#endif 



#define SAMPLEPERSECOND 44100 
#define WINDOWSIZE 1024 

// calculates the frequency resolution of audio signals 

#define DELTAFREQUENCY (double (44100) / double(1024)) 

// convert sample number to time in seconds. 

#define sample2time(s) (double(s) / doubl e(SAMPLEPERSECOND)) 



using namespace std; 



// 



// 
// 
// 
// 
// 
// 
// 
// 
// 
// 



signal — > o Delay — + Sum — + e 

+-> o 



I 



noise -+-> o + Adaptive 

I filter 
I /_.__ 

I 



— + 



I 

-+ 



const float STEP = O.Olf; 
const int NTAPS = 2048; 
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const int DELAY = NTAPS / 2; 



int 

cancel Noise(string waveFileName) 
int n; 

float *taps = nspsMallocC NTAPS ); 
float *dlyl = nspsMallocC NTAPS ); 
float err = 0; 



Readwave Fi 1 e wavFi 1 e (wave Fi 1 eName . c_s t r () ) ; 

int noOfSamplesPerChannel = wavFile.getNoOfSamplesPerChannel () ; 
Wri tewaveFi 1 e wavFi 1 eLef tRi ght ("h : /I ef t ri ght . wav" , 
wavFile.wf .nSamplesPerSec, true) ; 

cout « wavFile.wf .nSamplesPerSec « "\t M « noOfSamplesPerChannel 

« endl ; 



clock_t startTime = clock(); 

NSPLmsTapState tapStPtr; 
NSPLmsDlyState dlyStPtr; 

/// init taps delay line to zero 

for( n=0; n<NTAPS; ++n ) taps[n] = dlyl[n] = 0; 

/// initialize filter 

nspsLmsllnit( NSP_LmsDefaul t , taps, NTAPS, STEP, O.Of, 0, AtapStPtr 
/// Initialize delay line 

nspsLmsllnitDlyK AtapStPtr, dlyl, TRUE, &dlyStPtr ); 

/// Filter LEN samples using single-rate adaptive filtering 

// each pass process (one+delta) seconds of samples 
float *left = new float [noOfSamplesPerChannel ] ; 
float *right = new float [noOfSamplesPerChannel ] ; 
wavFile.getSamples(left , right, noOfSamplesPerChannel) ; 

for (int s=DELAY ; s<noOfSamplesPerChannel ; s++) { 

float d2 = nspsLmsl( AtapStPtr, AdlyStPtr, left [s]/32768.0f , 

float d = right [s-DELAY]/32768. Of ; 
err = d - d2; 

wavFi 1 eLef tRi ght. addsample(err*32768.f ,err*32768.f) ; 
//wavFi 1 eLef tRi qht . addsampl e (e r r , ri ght [s-DELAY] ) ; 
if (s%wavFile.wf .nSamplesPerSec == 0) 

cout « s/wavFile.wf .nSamplesPerSec « endl; 

} 



clock_t stopTime = clock(); 

cout « "Total processing time: " « (stopTime - 
startTime)/CLOCKS_PER_SEC « endl; 
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); 



err ) ; 



main 

nspFree( dlyl ); 
nspFree( taps ) ; 
delete[] left; deleteC] right; 

return 0; 

} 

// 



// - calculate total perceptual loudness ("gesamtLauthei t") 

// 



class Convol veAudioBuffer { 
public: 

static float* mask; 

static int masksize; 

static int noOfSamplesProcessed; 

public: 

Convol veAudioBuffer () {} 
-Convol veAudioBuffer () {} 



float operator () (float* windowSamples , int windowsize, int 
samplingRate, float samplescale) 

int windowSize2 = windowsize / 2; 
int windowSizeBase2 = 0; 
int temp = windowsize; 

while (temp > 1) { temp »= 1; wi ndowSizeBase2++; } 
float* result = new float [windowsize + masksize]; 
nspsConv(windowSamples, windowsize, mask, windowsize, 

result) ; 

float resultvalue = result [masksize-l] ; 
delete result; 

noOfSamplesProcessed++; 
return resultvalue; 

} 

}; 

float* Convol veAudioBuffer: :mask=0; 

int Convol veAudioBuffer: :masksize=1024; 

int convol veAudioBuffer: :noOfsamplesProcessed=0; 



// 



// - write audio test file to test mic's frequency response 

// 
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int 

writeTestwaveFile(int argc, char** argv) 

writeWaveFile wavFi le( M c: /temp/test .wav" , samplepersecond) ; 
double frequency = deltafrequency/2 ; // of the sinusoid 
double time = 0; 

int j=0; 

for (int loop=0 ; loop<WlNDOWSIZE/2 ; loop++, frequency+= 
DELTAFREQUENCY) { 

// cut of all frequencies above 15000 

if ((frequency-DELTAFREQUENCY/2) > 15000) break; 

cout « "Frequency: " « frequency « endl ; 

for (int i=0 ; i <samplepersecond*4 ; i++, j++) { 
time = SAMPLE2TlME(i) ; 
wavFile.addSample(short(10000 * 
sin(f requency*2*M_Pl*time)) ); 

for (int i=0 ; i <SAMPLEPERSECOND ; i++, j++) { 
wavFi 1 e . addsampl e (sho rt (0) ) ; 

} 

} 

return 0; 

} 

// 



// - get loudness of the left and right audio channel 

// 



bool 

cjenerateStatistics(int argc, char** argv) 

ReadwaveFile wavFile("d :/TestSignal ForMicrophone.wav") ; 
int s = wavFile.getNoOfSamplesPerChannel () ; 
short left, right; 

clock_t startTime = clock(); 

vpl : :AudioBuffer<vpl : :gesamtLauthei t> audiobufferl; 
vpl : :AudioBuffer<vpl : :qesamtLauthei t> audiobuffer2 ; 
for (int i=0 ; i<wavFile.getNoOfsamplesPerchannel () ; i++) { 
wavFi 1 e . qetsampl e (1 eft , ri ght) ; 
audiobufferl.addsample(left) ; 
audi obuffer2 .addsampl e(ri ght) ; 
if (i%wavFile.wf .nSamplesPerSec == 0) 

cout « i/wavFile.wf .nSamplesPerSec « "\t" « left 
« M \t" « right « endl; 

clock_t stopTime = clock(); 

cout « "Total processing time: " « (stopTime - 
startTime)/CLOCKS_PER_SEC « endl; 
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return 0; 

} 

// 



// - get loudness of the left and right audio channel 

// 



bool 

cjenerateStatistics2(int argc, char** argv) 

const int WINDOWSIZE2 = WINDOWSIZE / 2; 

ReadwaveFile wavFile("d:/Testsignal ForMicrophone.wav") ; 
int s = wavFile.getNoOfSamplesPerChannel () ; 

int windowSizeBase2 = 0; 

int temp = WINDOWSIZE; while (temp > 1) { temp »= 1; 
windowSizeBase2++; } 

float left[wiNDOWSiZE+2] , right [windowsize+2] ; 



// 



// CREATE CONVERSION MASK 

// 



const int MASKSIZE = 1024; 

const int MASKSIZE2 = MASKSIZE / 2; 

float maskD[MASKSIZE+2] ; nspsbset(0.0f , maskD, MASKSIZE+2) ; 
float maskMic[MASKSIZE+2] ; nspsbSet(0.0f , maskMic, MASKSIZE+2) ; 
// 



clock_t startTime = clockO; 

// read chunk at the beginning of wave file: 393*1024 samples = 
9*44100+5532 

for (int i=0 ; i<393 ; i++) wavFile .getsamples(left , right, 

WINDOWSIZE) ; 



FILE* testData = fopen("d:/mics. test" , "w") ; 

const int noOf Samples = 5 * 44100-183; // the length of one sinusoid 
int noOf Headi ngwi ndows = 44100/2 / WINDOWSIZE; // skip the first 
half second 

int noofwi ndows = 44100*3/WINDOWSIZE - noOf Headi ngwi ndows ; 

int noOfTail windows = noOfsamples/wiNDOWSlZE - noofwi ndows - 
noOf Headi ngwi ndows ; 

int noOfTailwi ndowsRest = noOfSamples - 
(noOf Headi ngwi ndows+noOf Wi ndows+noOfTai 1 wi ndows) *WINDOWSIZE ; 
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for (int i =0 , j =0 ; j<222 && i<wavFile,getNoOfSamplesPerChannel () 
i+=no0f Samples, j++) { 

cout « j « " : "; 

// process the current frequency 
complex<double> ratioSum = 0.0; 
int countNoOfRatios = 0; 

for (int x=0 ; x<noOfHeadingwindows ; x++) 
wavFi 1 e . getsampl es (1 eft , ri ght , WINDOWSIZE) ; 

for (int n=0 ; n<noOfwindows ; n++) { 

wavFi 1 e . getsampl es (1 ef t , ri ght , WINDOWSIZE) ; 

// calculate frequency magnitudes 
nspswi nHammi ng (1 eft , WINDOWSIZE) ; 
nspsRealFft(left, 10, NSP_Forw) ; 
nspswi nHammi ng(ri ght , WINDOWSIZE) ; 
nspsRealFft(nght , 10, NSP_Forw); 

float windowMagnitudeL[wiNDOWSlZE2+l] ; 
nspcbMag((const SCplx*) left, windowMagni tudeL, 

float windowMagni tudeR[WlNDOWSlZE2+l] ; 
nspcbMag((const SCplx*) right, windowMagni tudeR, 



WINDOWSIZE2+1) ; 
WINDOWSIZE2+1) ; 



// CREATE MICROPHONE STATISTIC 

float* pmaxL = std: :max_el ement (windowMagni tudeL, 
windowMagni tudeL+WINDOWSlZE2) ; 

float* pmaxR = std: :max_el ement (windowMagni tudeR, 
windowMagni tudeR+wiNDOWSlZE2) ; 

int imaxL = int(pmaxL - wi ndowMagni tudeL) ; 

int imaxR = int(pmaxR - wi ndowMagni tudeR) ; 

//float minv = *std: : mi n_el ement (wi ndowMagni tude, 
wi ndowMagni tude+wi ndowsi ze2) ; 

//fprintf(testData, M %d %d %f %f %f %f %f %f\n M , 

imaxL, imaxR, 

// *pmaxL, *pmaxR, 

// left[imaxL*2] , left[imaxL*2+l] , 

right[imaxR*2] , right [imaxR*2+l]) ; 

complex<float> r(right [imaxR*2] , right [imaxR*2+l]) ; 
complex<float> 1 (left[imaxL*2] , left [imaxL*2+l]) ; 
complex<float> ratio = 1/r; 
if ((imaxL == j) && (imaxR == j)) { 

ratioSum += ratio; 

countNoof Rati os++ ; 

//fprintfCtestoata, "%d %d %f %f %f %f\n", 
imaxL, imaxR, *pmaxL, *pmaxR, ratio. real () , ratio. i mag ()) ; 

} 

} 

ratioSum /= countNoOfRatios; 

fprintf (testData, "%f %f\n", ratioSum. real () , 
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rati oSum.i mag ()) ; 

maskD[j*2] = rati oSum. real () ; 
maskD[]*2+l] = rati oSum.i mag O ; 

for (int x=0 ; x<noOfTailwi ndows ; x++) 
wavFile.getSamples(left, right, WINDOWSIZE) ; 

wavFi 1 e . getSampl es (1 eft , ri ght , noOfTai 1 wi ndowsRest) ; 

} 

fclose(testData) ; 

clock_t stopTime = clockO; 

cout « "Total processing time: " « (stopTime - 
startTime)/CL0CKS_PER_SEC « end!; 



// 



// CREATE CONVERSION MASK 

// 



// the first 3 frequencies are unsave 

for (int x=0; x<3 ; x++) {maskD[2*x]=l; maskD[2*x+l]=0; } 
for (int x=222 ; x<MASKSl2E2 ; x++) maskD[2*x]=l; 
of stream outPut2("conversionMaskOrg. txt") ; 

copy(maskD, maskD+MASKSiZE, ostrearrui terator<float>(outPut2 , "\n")); 
outPut2 ,close() ; 



nspsCcsFft(maskD, 10, NSP_lnv) ; 
ofstream outPut( M conversionMask. txt") ; 

for (int i=0 ; i<MASKSlZE ; i++) { maskMic[(i+MASKSlZE2-l) % 
MASKSIZE] = maskD[i]; } 

maskMic[MASKSlZE-l] = 0; 

copy(maskMi c, maskMic+MASKSlZE-1, ostream_i terator<float>(outPut , 

"\n")); 

// determine the cut-off coefficients 
/* float* posln = find_if (maskMic, maskMic+MASKSlZE, 

bind2nd(greater<float>() ,0.01f)) ; 

int usedMaskSizel - posln - maskMic; 

int usedMasksize3 = MASKSIZE-1 - usedMaskSizel -1; 

for (int i=0 ; i<MASKSlZE-usedMaskSizel ; i++) { maskMic [i] = 
maskMic [i+usedMaskSizel] ;} 

int usedMaskSize = usedMaskSize3 - usedMaskSizel + 1; 

copy (maskMi c , maskMi c+usedMaskSi ze , ost ream_i terator<f 1 oat>(outPut , 

"\n")); 

*/ outPut-flush() ; 
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return 0; 



// 



// - correlate the response of the 2 mics to the sinusoid of a certain 
frequency 

// 



bool 

correlateChannelsO 

ReadwaveFile wavFile("d :/TestSignal ForMicrophone .wav") ; 
int s = wavFile.getNoOfSamplesPerChannel () ; 

const int noof Samples = 5 * 44100-183; // the length of one sinusoid 
float* left = new float [noof Samples] ; 
float* right = new float [noof samples] ; 

clock_t startTime = clock(); 

// read chunk at the beginning of wave file: 393*1024 samples = 
9*44100+5532 

wavFi le. getsampl es (left, right, 3 * 44100); 

wavFi le. getsampl es (left , right, 3 * 44100); 

wavFi le. getsampl es (left, right, 3 * 44100); 

wavFi le. getsampl es (left, right, 5532) ; 

wavFi 1 e . getsampl es (1 eft , ri ght , noof Sampl es) ; 

for (int i=0,j=l ; j<222 && i<wavFi le.getNoOf Sampl esPerchannel () ; 
i+=noOf Samples, j++) { 

cout « j « ": "; 

// REMEMBER: cut of first half second and last half second 
// find best correlation 

wavFi 1 e . getsampl es (1 eft , ri ght , noof Sampl es) ; 
char name [1000] ; sprintf (name, n d:/%03d . wav" , j); 
WritewaveFile wavFreq(name, SAMPLEPERSECOND, true); 
wavFreq . addsampl es (1 eft , ri ght , noof Sampl es) ; 
//avFreq.addSamples(left+44100/2, right+44100/2 , 2*44100); 

for (int n=0 ; n<noOf Samples ; n++) left[n] = -left[n]; 
const float stepFrequency = float(44100) / float(1024); 
int noofsamplesperwave = int (44100 / (j * stepFrequency + 
stepFrequency/2 .Of)) ; 

cout « noofsamplesperwave « "\t"; 

float *result = new float[201]; // [noof Sampl esPe rwave+l] ; 
//nspsCrossCorr(left+44100/2, 3*44100, right+44100/2, 
3*44100, result, 0, noofsamplesperwave); 

nspsCrossCorr(left+44100/2, 2*44100, right+44100/2, 2*44100, 

result, 0, 200); 

//float* minPos = min_element(result , 
resul t+noOf Sampl esPerwave+1) ; 

// return local maxs 

for (int n=l ; n<200-l ; n++) { 

if ((resul t[n-l] <= result[n]) && (result[n+l] <= 
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result[n])) 

cout « n « "\t"; 

} 

delete result; 
cout « end! ; 

clock_t stopTime = clock(); 

cout « "Total processing time: " « (stopTime - 
startTime)/CLOCKS_PER_SEC « endl; 

return 0; 

} 



// 



// - diente dazu, Probleme mit der FFT auszuraeumen 

// 



void 

cjenerateFFTTestPatterns(int argc, char** argv) 

const int MASKSIZE = 1024; 

const int MASKSIZE2 = MASKSIZE / 2; 

float maskD [MASKSIZE+2] ; nspsbSet(0.0f , maskD, masksize+2) ; 

of st ream outPutorg ("conversionMaskOrg.txt") ; 
nspsbset(0.0f , maskD, MASKSIZE+2) ; 
for (int i=0 ; i<MASKSlZE2 ; i++) { 

float iFloat = float(i) / float (MASKSIZE2); 

maskD [i] = exp(-iFloat*i Float / (2 * 0.0001)); 

maskD [MASKSiZE-l-i] = exp(-i Float*i Float / (2 * 0.0001)); 

} 

for (int i=0 ; i<MASKSlZE ; i++) { 

//outPutorg « maskD [i] « endl; 

} 



nspsReal Fft(maskD, 10, NSP_Forw) ; 

//float mag [MASKSIZE2+2] ; nspsbSet(0.0f , mag, MASKSIZE2+2) ; 
//nspcbMag((const SCplx*)maskD, mag, 512); 
for (int i=0 ; i<MASKSlZE ; i++) { 

outPutorg « maskD [i] « endl; 

//outPutorg « mag[i] « endl; 

} } 

// 



// - diente dazu, Probleme mit der fft auszuraeumen 

// 
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void 

cjenerateFFTTestPatterns2(int argc, char** argv) 

const int MASKSIZE = 1024; 

const int MASKSIZE2 = MASKSIZE / 2; 

float maskD[MASKSlZE+2] ; nspsbSet(0.0f , maskD, MASKSIZE+2) ; 
for (int i=0 ; i<MASKSiZE2 ; i++) maskD[2*i] = l.Of; 



of stream outPutOrg("conversionMaskOrg. txt") ; 
nspsCcsFft (maskD, 10, NSP_lnv) ; 
for (int i=0 ; i<MASKSlZE ; i++) 

outPutOrg « maskD [i] « end!; 

} 



// 



// - Online Training 



// 



void 

onl i neTrai ni ng(st ri ng waveFi 1 eName) 

ReadWaveFile wavFile (waveFi 1 eName. c_str()) ; 

int noOfSamplesPerchannel = wavFile.getNoofsamplesPerChannel () ; 



Wri teWaveFi 1 e wavFi 1 el_ef tRi ght (" c : /temp/1 eft ri ght . wav" , 
SAMPLEPERSECOND, true); 

clock_t startTime = clock(); 

const int NOOFSAMPLES = 1024; 

const int NOOFSAMPLES2= NOOFSAMPLES/2 ; 

const int POWEROFSAMPLES = 10; 

const float threshold = 100 * NOOFSAMPLES; 

float left [NOOFSAMPLES+2] ; float leftMag [NOOFSAMPLES2+1] ; 
float right [NOOFSAMPLES+2]; float rightMag[NOOFSAMPLES2+l] ; 
float ratio [NOOFSAMPLES+2] ; nspsbzero(ratio, NOOFSAMPLES+2); 
int count [NOOFSAMPLES2+1] ; memset(count , 0, sizeof(int) * 

(NOOFSAMPLES2+1)) ; 

// each pass process (one+delta) seconds of samples 
for (int s=0 ; s<noofSamplesPerChannel /NOOFSAMPLES ; s++) { 
wavFile.getsamples(left, right, NOOFSAMPLES); 

nspsRealFft(left, POWEROFSAMPLES, NSP_Forw); 
nspsRealFft(right, POWEROFSAMPLES, NSP_Forw); 
nspcbMag((const SCplx *)left , leftMag, noofsamples/2+1) ; 
nspcbMag((const SCplx *)right, rightMag, NOOFSAMPLES/2+1) ; 

float leftMean = nspsMean(leftMag, NOOFSAMPLES2+1) ; 
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float rightMean = nspsMean(rightMag, NOOFSAMPLES2+1) ; 
for (int i=0 ; i<N00FSAMPLES2+l ; i++) { 

if C(leftMag[i] > threshold) && (rightMag[i] > 

threshold)) { 

// cout « i « " " ; 
count[i]++; 

complex<float> leftc (left [2*i], left 

[2*i+l]); 

compl ex<f 1 oat> ri ghtc(ri ght [2*i ] , 

right[2*i+l]) ; 

compl ex<float> ratioC = leftc / rightC; 
ratio[2*i] += rati oc. real () ; 
ratio[2*i+l] += ratioC. i mag () ; 

} 

} 

COUt « (s*NOOFSAMPLES)/SAMPLEPERSECOND « "\t" « leftMean 

« "\t" « rightMean « end! ; 



// Calculate deconvolution mask 
for (int i=0 ; i<NOOFSAMPLES2+l ; i++) { 
if (count [i]) { 

ratio[2*i] /= count [i]; 

ratio[2*i+l] /= count[i]; 

} else { 

ratio[i] =0.5; 
ratio[2*i+l] = 0; 

} 

} 



nspsCcsFft(ratio, POWEROFSAMPLES, NSP_Inv); 



ofstream outPut("conversionMask. txt") ; 
float maskMic [noofsamples] ; 

for (int i=0 ; i <NOOFSAMPLES ; i++) { maskMic [(i+NOOFSAMPLES2-l) % 
noofsamples] = ratio[i]; } 

copy(maskMi c , maskMi c+NOOFSAMPLES , ostream_i terator<f 1 oat> (output , 

"\n")); 

output .flush() ; 



clock_t stopTime = clock(); 

cout « "Total processing time: " « (stopTime - 
startTime)/CLOCKS_PER_SEC « endl; 
} 



// 



// - MAIN 



// 
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int 

main(int argc, char** argv) 

cancel Noi se("h :/openHouse . wav") ; 
//generateStati sties (argc, argv) ; 
return 0; 

//string wavlnputFileName = M d: /shortwave. wav" ; 
string wavlnputFileName = "d:/cap2 .wav" ; 
onl i neTrai ni ng(wavlnputFi 1 eName) ; 

// 



// CREATE CONVERSION MASK 

// 



const int MASKSIZE = 1024; 

const int MASKSIZE2 = MASKSIZE / 2; 

float maskD[MASKSlZE+2] ; nspsbset(0.0f , maskD, MASKSIZE+2); 
float maskMi c[MASKSIZE+2] ; nspsbSet(0.0f , maskMic, MASKSIZE+2); 
/* 

i f stream i nPut ("conversi onMaskorg . txt") ; 
assert(inPut != 0) ; 
for (int i=0 ; i<MASKSlZE ; i++) { 
if (! input. eof()) 

input » maskD[i] ; 

} 



nspsCcsFft(maskD, 10, NSP_Inv) ; 
ofstream outPut("conversionMask. txt") ; 

for (int i=0 ; i<MASKSlZE ; i++) { maskMi c[(i+MASKSlZE2-l) % 
MASKSIZE] = maskD [i]; } 

maskMic [MASKSIZE-1] = 0; 

V 

ifstream inPut("conversionMask.txt") ; 

copy(i stream_i terator<float>(inPut) , istream_i terator<float>() , 
maskMic) ; 

// determine the cut-off coefficients 
/* 

float* posln = find_if (maskMic, maskMi c+masksize, 
bind2nd(greater<float>() ,0.01f)) ; 

int usedMaskSizel = posln - maskMic; 

int usedMaskSize3 = MASKSIZE-1 - usedMaskSizel -1; 

for (int i=0 ; i<MASKSlZE-usedMaskSizel ; i++) { maskMic [i] = 
maskMi c [i +usedMaskSi zel] ; } 

int usedMasksize = usedMasksize3 - usedMaskSizel + 1; 

copy (maskMi c , maskMi c+usedMaskSi ze , ostream_i terator<f 1 oat>(outPut , 

"\n")); 

output. flush() ;*/ 

for (int i=0 ; i<1024/*usedMasksize*/ ; i++) { 
cout « maskMi c[i] « endl ; 

} 
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int usedMaskSize = 1024; 

Convol veAudi oBuf f e r : : mask=maskMi c ; 

convol veAudi oBuf f e r : : masksi ze=usedMasksi ze ; 



ReadWaveFi 1 e wavFi 1 e (wavlnputFi 1 eName . c_st r() ) ; 
//ReadwaveFi 1 e wavFi 1 e("d : /TestSi gnal ForMi crophone . wav") ; 
int s = wavFile.getNoOfsamplesPerchannel () ; 

int noofSamplesPerchannel = wavFile.getNoofSamplesPerchannel () ; 
wri teWaveFi I e wavFi 1 ei_ef tRi ght ("c : /temp/1 ef t ri ght . wav" , 
SAMPLEPERSECOND, true); 

clock_t startTime = clockQ; 



float left 
float right 
float result 



SAMPLEPERSECOND + MASKSIZE] ; 
SAMPLEPERSECOND + MASKSIZE] ; 
SAMPLEPERSECOND + 2 * MASKSIZE] ; 



// prefill buffers 

wavFi 1 e . getsampl es (1 eft , ri ght , MASKSIZE) ; 

// each pass process (one+delta) seconds of samples 
for (int s=0 ; s<10 ; s++) { 

wavFile-getSamples(left + MASKSIZE, right + MASKSIZE, 
SAMPLEPERSECOND) ; 

nspsConv(right, SAMPLEPERSECOND+MASKSIZE, maskMic, MASKSIZE, 

resul t) ; 

//nspsbMpyl(1.2f , right+MASKSIZE/2 , SAMPLEPERSECOND+1) ; 

//nspsbSut>2 (1 ef t+MASKSlZE/2 , resul t+MASKSlZE-1 , 
SAMPLEPERSECOND+1) ; 

wavFi 1 eLef tRi ght . addsampl es (1 ef t+MASKSlZE/2 , 
resul t+MASKSIZE-1, SAMPLEPERSECOND+1) ; 

//nspsbMpyl(1.0f/1.2f , right+MASKSIZE/2 , 
SAMPLEPERSECOND+1) ; 

memcpy(right , right+SAMPLEPERSECOND, sizeof (float) * 

MASKSIZE) ; 

memcpy(left, left +SAMPLEPERSECOND, sizeof (float) * 

MASKSIZE) ; 

cout « s « endl ; 

} 

clock_t stopTime = clock(); 

cout « "Total processing time: " « (stopTime - 
startTime)/CLOCKS_PER_SEC « endl; 



return 0; 
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